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ABSTRACT 


Ion acoustic phenomena are studied in a cylindrical geometry 
for two distinct cases. A large amplitude compressive pulse is 
seen to evolve into solitons. The evolution of these solitons 
and their dependence on initial conditions show a similarity to 
previous work on one dimensional solitons. Dimensionless scaling 
arguments are used to distinguish between the two cases . In the 
presence of a steady state uniform cylindrical beam, approximated 
by a ring in an ion-ion beam instability is observed. This 

instability exists- for a limited range of beam velocities and shovrs 
a marked similarity to the strictly one dimensional ion-ion beam 
instability. Solution of the appropriate dispersion relation 
shows agreement with the observed phenomenon. 
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I. INTRODUCTION 

The investigation of ion acoiostic waves has been the sub- 
ject of extensive study in recent years. With the invention of 
the double plasma (DP) device at UCLA [Taylor I969A; Taylor 1972] 
new and wider areas of investigation were opened. Large ampli- 
tude planar waves can be efficiently launched in a plasma with 
variable electron to ion temperature ratio [Mackenzie 19Td]» 

Steady state large diameter ion beams can also be injected into 
the plasma. 

The initial work with a DP device was an investigation of 
collisionless electrostatic shocks [Taylor 1969B]. Ramp 
excitations of variable height and rise time were applied to the 
plasma. For large amplitude excitations ( “ < ,5) the initial 
ramp was observed to steepen and form a collisionless electro- 
static shock. A small fraction of ions were reflected by the 
potential Jump of the shock forming a double humped ion distri- 
bution in front of the shock. 

These reflected ions represent a localized ion beam in the 
plasma. Stem and Decker [1971] report observation of an instability 
due to these reflected ions. A velocity gradient was established 
by injecting an ion beam from the target into the driver chamber. 
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When a large ait^^litude pulse is launched into the target chamber 
the reflected ions travel balistically from the point of reflection 
until the velocity gradient brings the precursors into resonance 
with the local ion acoustic velocity. A resonant transfer of 
energy from precursors to noise fluctuations results in the growth 
of the noise until it engulfs and destroys the original wave. 

A similar study of turbulence associated with precursors 
of ion acoustic shocks by Means et al- [1975] attributes the observed 
turbulence, in their case, to the three dimensional nature of the 
instability. With no velocity gradient present there are modes 
oblique to the direction of the wave front which resonate with the 
reflected precursors . It is these off axis modes that result in 
the turbulence. Three qualitatively different effects were observed. 
For high electron to ion temperature ratios (T^./T; ) and small ampli- 
tude, laminar shocks with no tiorbulence were observed. For moderate 
Te/Ti ratios or amplitude turbulence, restricted to the region occupied 

by the precursor ions was observed. For low T /t.. or large ampli- 

e 1 

tude, the turbulence engulfed the shock and destroyed its structure. 

An investigation of the heating of an ion beam was made by 
Taylor and Coronitti [1972]. A uniform planar ion beam with velo- . 
city greater than twice the ion acoustic velocity was injected 
from the driver chamber into the target chamber. The ion- ion 
beam instability excited off-axis inodes that propagated oblique 
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to the beam direction. The ion beam is seen to be heated ijntil it 
forms a quasilinear plateau on the background ion distribution. 

The evolution of non-linear compressive pulses that are governed 
by the Korteweg-de Vries equation -was made possible by the DP 
device. Ikezi, Taylor and Baker [1970] investigated the properties 
of single solitons . The width, amplitude and velocity were shown 
to be consistent with the definition of solitons. The interaction 
of two solitons was also studied. Hershkowitz, Romesser and 
Montgomery [1972] compared the evolution of an initial condition 
into multiple solitons with the analytic theory. A s5.milar report 
has been made by Ikezi [ 1973 ] in which both pulse and continuous 
wave initial conditions are treated. 

The objective of the experiments reported in this disser- 
tation was to investigate in a cylindrical geometry two kinds of 
ion acoustic phenomena. A study of the evolution of weakly non- 
linear compressive pulses in a cylindrical geometry was made to 
determine if solitons or soliton-like objects existed. Once their 
existence was determined a detailed comparison between planar and 
cylindrical solitons became possible. The second area of investi- 
gation is the ion-ion beam instability in a cylindrical geometry. 

For steady state radially ingoing beams an ion-ion beam instability 
is expected to occur. The dependence of the instability on beam 
velocity and beam density has been investigated and compared with 
the analytic theory for a linear instability. 
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II. EXPERIMENTAL APPARATUS 

The experiments reported here were performed in the University 
of Iowa double plasma device. A modification was made to allow oper- 
ation in a cylindrical geometry. The University of Iowa double 
plasma device (Fig. 1) consists of two aluminum cylinders (length 
37 CKi; diameter 38 cm) placed end to end but separated by an insu- 
lating ring. In each chamber 36 filaments are located symmetrically 
near the outside wall. Two cylindrical fine mesh screen (200 lines/ 
inch) of radii 10 and 11 cm and length ~ 30 cm are mounted concen- 
trically in one end of the device. Electrons emitted from the filaments 
which are biased at -4o volts with respect to the wall produce the plasma. 
The outer cylindrical screen is allowed to float causing it to charge 
up negatively preventing the flow of electrons between the inner and 
outer plasma. The inner cylindrical screen is grounded to help de- 
termine the plasma potential of the interior plasma. A shield was 
used to prevent line of sight paths between the outer (driver) plasma 

and the inner (target) plasma. Typical operating parameters are 
/ \ 9-3 

ion density ^ 10 cm , electron temperature (T^) 1-5.0 eV, 

and an ion temperature (T^) ^ .2 ev [Hudis et al . I 968 ] at a pressure 
, -4 

of 3 X 10 Torr. 
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Cooling coils are wrapped around the chamber and attached 
using Stycast 285 0 FT epcocy. The system is evacuated using a 4 inch 
oil diffusion pump to an ultimate base pressure of ~ 5 X 10 torr 
after pumping down. This assures us of negligible impurity con- 
centration if data is not taJcen until the ultimate base pressure 
is reached. 

There are two modes of operation for the DP device. If the 
potential of the driver chamber wall is raised with respect to the 
target chamber v/all (ground) a uniform symmetric cylindrical beam 
is launched. The beam energy is adjustable^ being a function of , the 
applied potential difference. The beam is detected using a gridded 
energy analyzer as shown in Fig. 2. The front screen (200 lines/inch) 
is allowed to float and prevents electrons from entering the analyzer. 
The collector plate is biased at -67.5 volts to collect ion current. 
The voltage to the center grid is ramped with a sawtooth voltage. 

Only those particles with sufficient directed momentum to overcome 
the barrier represented by the second grid can be collected. This 
particular design gives directional resolution as will be seen later. 
From geometrical considerations the analyzer can detect particles 
in a cone of half angle Id . The fields generated in the plasma 
wiUL tend to broaden this. An energy resolution of .1-.2 eV is 
possible but deterioration of resolution occurs with use and exposure 
to the plasma. Figure 5 shows typical energy analyzer traces for 

¥r 

Emerson and Cumings, Inc. Canton, Mass, 
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several beam energies. The signal represents the flux of particles 

"with energy greater than ^ plotted versus f . If we then differen- 

s s 

tiate with respect to time (assuming $ is linear with time) we 

s 

recover the velocity distribution function f (v) plotted versus . 

Figure 4 shows a differentiated energy analyzer trace. The temper- 
ature of the background ions is seen to be ajjproximately .2 eV and 
a well defined beam is seen with energy 3.2 eV. The width of the 
beam appears slightly wider than that of the background ions . To 
determine the temperature, however, we must replot the data as f(v) versus 
V. Figure 5 shows the effect of doing this. The beam as expected 
will be much colder after being accelerated. To find the beam 
density we must look at the data that has been replotted as f(v) vs v. 
Measurement of electron density is made using a Langmuir 
probe. To determine the density a large cylindrical probe (length 
1 inch; diajiieter l/8 inch) is used. For detection of waves a small 
wire biased above the plasma potential is used to detect electron 
satxiration current. 

Before looking at some typical data it must be noted that it 
is virtually impossible to form truly cylindrical screens. There is 
most probably a variation in the center of curvature of ~ .5 cm in 
the screen and for this reason no divergence is observed at r = 0. 

In addition there is a competition between damping and geometrical growth. 
Experimentally, they are found to be very closely balanced. 

Upon application of a half sine wave pulse to the driver 
chamber a eylindrically symmetric ingoing wave is launched at r = 10 cm. 
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The -wave propagates inward ^ converges at the center and emerges as 
a cylindrically synmetric outgoing wave. If we set the probe at a 
fixed position and loolc at the perturbed electron density as a function 
of time, we see first the ingoing wave and at a later time the wave 
after it has converged at the center and emerged as an outgoing 
wave . This is apparent in Figure 6. where the perturbed electron 
density as a function of time is shown for 1 cm intervals. The ion 

c 

acoustic velocity as measured from this is 2.l6 X 10"^ cm/sec which 
indicates an electron temperature of ~ 2.0 eV. 

The presence of both ingoing and outgoing waves makes it 
difficult to make a detailed comparison between cylindrical and planar 
phenomenon. In order to make this comparison a major modification 
was made. The cylindrical screens were removed and a rectangular 
chamber was added to the end plate as in Fig. ^ deformable sepa- 
ration screen that was adjustable from a plane to a half cylindei” of 
radius 10 cm was used. Then, under nearly identical plasma conditions, 

a comparison of planar waves and half cylindrical waves was made. 

* 

Using a lock-in amplifier it is possible to take an inter- 
ferometer trace of a continuous wave excitation. Figure 8 shows 
the results for a planar separation screen with a frequency of 100 kHz. 
The lock-in amplifier detects the amplitude of the wave at a given 
phase with respect to the applied signal. For a plane wave we expect 
it to have the form: 


Par model 


HR -8. 
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i(k + ik. )r + 

A(r)=ReAe ^ ^ 

o 


-k r i(k r + p) 

= A e ^ Re e 
o 


-k.r 

= A e ^ cos(k r + 0) 
o \ p r/ 


For the case shovm in Figure 8 we can measure k and k. . From this 
k. r 1 

we find ^ = .o44. From this and assuming ion Landau dumping we 
^r 

can determine T /t. to be approximately 15. This agrees with T ~ 1 eV 

^ dL 6 

and T . .1 eV, 

1 

A similar interferometer trace for the case of the separation 
screen deformed to a half cylinder of radius 10 cm is shown in 
Figure 9» Unlike the planar case, the wavelength is no longer 
constant for the cylindrical waves . The received amplitude for 
the cylindrical waves will be the sum of J^(kr) and Y^(kr), 


A(r) == Re{A^[J^(kr + 0) + iY^(kr + 


where k = k + ik. . For the case 0 = 0 k can be expressed as 
r 1 ^ r ^ 

|k|re^ . Tables of this are available [Nat. Comp. Lab., 19^7^ 1950], 



Figure 10 shows a conputed interferometer trace for cylindrical 
waves using the values and as determined from the planar 
interferometer traces. The wave length is not constant in this 
case^ though the disparity is much less than is measured from 
Fig. 9. 

Data is taken using a Biomat ion 6lOB transient recorder 
which acts as a fast analog to digital converter, A description 
of this system is given in the Appendix. 
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III. SOLITOl^S 
A . Theory 

The properties of systems that support soliton solutions has 
generated a great deal of interest in recent years [Scott et al . 1975] 
Extensive work in plasma physics has been carried out since Washimi 
and Taniuti [19^6] showed that weakly non-linear ion acoustic waves 
are governed by the Korteweg-de Vries (KdV) equation; 


ht dx 2^5 


( 1 ) 


Equation 1 is written in a frame moving at the ion acoustic velocity. 


n is the fractional density perturbation, x is the distance in the 

r kT I 1/2 

moving frame measured in electron Debye lengths = — k ^ 

Umi e^l 

t- ♦ o 


and 


t is measured in ion plasma periods 


■i- 

^ \J^Tm 


2 * 


The asymptotic stationary solution corresponds to a super- 
position of spatially separated solitary pulses, solitons. 


n(x, t) = y. 2e ^ Sech' 




c is the ion acoustic velocity , e. and 0. are constants 



m. J 'J 
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determined by the solution of an appropriate time independent 

Schrodinger equation with a potential well that is proportional 

to -1 times the initial spatial density perturbation [Hershkowitz 

et al . 1972 ]♦ The as^ptotic stationary soliton solution arises from 

a competition between non-linearity and dispersion. 

There are several properties of one dimensional solitons 

that can be used to identify them. Among them are; l) Arbitrary 

positive (compressive) density perturbations evolve after sufficient 

time into a superposition of spatially separated solitons, or solitary 

pulses; 2) The number and amplitude of the solitons is determined 

by the solution of an appropriate time independent Schrodinger 

equation: 5) The soliton velocity is given by [1 + ~ (”) ]c 

5 n max s 

where (— ) is the maximimi density perturbation of the soliton: 
n max 

4) The spatial width of the soliton. is proportional to [(^) ] 

' n max 

which implies that the product of the square of the maximum ampli- 
tude multiplied by the width is a constant; 5) Solitons retain their 
identity upon collision with other solitons. Experimentally these 
properties have been verified by Ikezi ^ [1970] and Hershkowitz 

^ [1972]. 

Maxon and Vlecelli [197^A] have reported the first work 
on solitons of dimensionality greater than one. They considered 
spherically symmetric perturbations and derived a modified KdV 
equation appropriate to a spherical geometry: 
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t ^ 


Sn , 1 iPn 

5 ? 2 ^ 5 = 


0 


This equation is valid except in a region near the origin. A simi- 
lar equation can be derived for cylindrical coordinates [Maxon and 
Viecelli 197^ B] starting with the fluid equations appropriate for 

a cylindrical geometry with no z or dependence. 


+ ::;^(nv) + — 0 Continuity equation (2) 

ot or r 


^v hv . hf 
ht ^ hr ~ ^r 


Momentum equation 


O) 


2 r hr 
^r 


n 


Poisson’s equation 


ih) 


kT 

In these equations $ is in units of — - — , n in units of n , r in units 

6 O 

of , t in units of u) ^ and v is in units of c . 

De pi s 

We now make the transformation to a frame moving inward at 
the ion acoustic velocity. 


§ = -e^/^(r + t) 



15 



where e is a small parameter that will he \jsed for expansion. The 
scaling used in the transformation is chosen to be the same as 
that of a soliton. We now expand the potential density and velocity 
in powers of e. 


•, ( 1 ) 2 ( 2 ) 

n = l+en +en^ ^ 

V = 

§ ~ 

The time and space derivatives transform to 

± - il A 4. 2? A - A/2 a ^ .5/2 a 

St St S| 5T! Stl S? ST) 

A=A.A + ®"A - eVs A. 

Sr Sr S5 Sir 6T! ~ ■ SJ 


Prom the transformation we find that 



r.= -e 

Equations 2, 3 
respectively; 

^5/2 A 
ST) 


^5/2 A 
ST) 


Xh 


+ T)) 


and 4 then transform to the following equations 






+ en 




p+ni ^ ^ 

(5) 

[vW . cv(2)] . .3/2 S[,(l) , ,^(2)3 

5s 

- e5/V"^ ^ 

= ^ ( 6 ) 



2 ii r il) , ,j(2), A 

e ^ L® + e« ] + (£g + n) 31 Ll + ef ] 


,,(1) 2,(2) , (1) ,2 (2) 

= e$' ' + G $' ^ + 7T €n' -€n 


(?) 


We can now collect terms of various orders of e. 
0(e) from Eq. 6. 




0(e^/^) 

65 6g 


for 0 as |^j -> 
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Then we obtain 

= 0 

yM + ^ 0 

To find the differential equation governing n^^^, and v^^^ we 

must go to higher order. 

O(e^) 


^ $( 2 ) ^ . „( 2 ) 

3|^ ^ 

To this point we have not had to account for contributions 
1 1 

arising from terms- with — dependence. The — factor gives rise to 

r r 

a factor of the form 



JlL 

eg ■*■ T1 

The order of this term depends upon the denominator. If <k j T|| 





IT 


this is of while if |ti| « |ef| this is of The 

former will be assumed and its implications will be discussed later. 


0(e5/2) 


3n^^^ 3 r_(l) (l)-, 

~a^ ^ ■ a? ^ 


il[v<2) 

35 



0 


( 9 ) 


3v( 2) av(^) 3$(2) 


( 10 ) 


( 2 ) 

Solving Eg. 8 for n' ' anii Eg. 10 for — rr — and substituting these 
results into Eg. 9 obtain: 


. .( 1 ) 1 „ 

which is the modified KdV eguation appropriate for cylindrically 
symmetric ingoing perturbation. 

In deriving this it was assumed that je^l « |Tj|. To under- 
stand the implications of this we must consider how the initial 
condition is entered into the problem. The initial condition is 
centered at ^ = 0 in the moving frame. 



iB 


= -€^^^(r + t ) = 


0 


r = -t 
o o 




where r is the radius at which the initial perturbation is applied, 
o 

The initial condition is then specified as initial 

time is negative and time increases as the pu3.se propagates inward. 
The time T] = 0 at which je^j « \'\\\ cannot possibly hold is the time 
for an ion acoustic wave to propagate to the origin where an expan- 
sion in orders of e would not be valid in any case . 

The value of g represents the distance traveled in the 
moving frame. The solitoiis have a velocity that is dependent upon 
the amplitude and in addition we expect an approximately growth. 


Assuming this we can approximately describe the soliton velocity as: 


A 6n 
V = A --- 

n o 


( "o V/2 


r - c t‘ 
o s 


where t^' = 0. where r = r and — 1 . is the soliton amplitude at 
^ o n o 



r = r^. The constant A describes the detailed velocity amplitude 
dependence and is of order 1. The distance in the moving frame is 
then given by 




The quantity T} is given by 


R 



The relation 


hsi « hi 


then becomes 


ML 

n o 


,r - c t' 
\ o s 


1/2 R 

s c 


We can now choose representative values for t' and evaluate the 

r 

inequality. For t ^ ^ which represents the time for an ingoing 

s 

linear wave to propagate half the distance to the origin we obtain; 
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A ^1 « 1 

n ' o 


which is true for moderate —I . 
3r 

For t^ = 

s 

we obtain 


A — I « I 
n 'o D 


Because of the presence of damping the criterion is less stringent, 

since a — ^ growth is not observed. 

- r 

The evolution of a weakly non-linear compressive pulse in 
the cylindrically symmetric geometry should be governed by the 
derived equation. This equation differs from the equation appro- 
priate to spherical geometry in the factor of l/2 associated with 





B. Presentation of Measurements 
The propagation of low amplitude pulses in the fully sym- 
metric geometry has already been noted. For low amplitude pulses 
the ingoing wave and outgoing wave are similar in shape and width 
to the applied pulse. As the amplitude is increased this is no 
longer true. Figure 11 shows the received signal detected at 
r = 0.5 cm as the height of the applied pulse is increased. For low 
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amplitudes a linear ion acoustic wave is launched. As the amplitude 

is increased the wave grows in amplitude, steepens and forms into an 

object that can be identified as a soliton [Hershkowitz and Romesser 

1974]. For still higher applied voltages the soliton decreases in artipli 

tude and a new faster pulse appears. This new pulse is a burst of 

particles, or a pseudowave [Alexeff et al . I968], traveling balisti- 

cally through the plasma with energy equal to the applied potential 

kT 

pulse. For pulse amplitudes less than ~ ^ no pseudowave is launched. 

e 

Figure 12 shows a comparison of the propagation of a linear 

< Vjo) ion acoustic wave and a large amplitude non-linear 
n 

o. Vflo) wave. At r = 9 cm for the linear pulse we see first in 

time an ingoing wave that closely resembles the applied half sine 

wave pulse, - followed at a later time by a similar outgoing pulse. 

As the probe is moved closer to the center the ingoing and outgoing 

pulse appear closer in time and finally merge at r = ,5 cm. At this 

point the received signal again closely resembles the applied half 

sine wave pulse. For the laxge amplitude pulse at r = 9 cm the 

ingoing wave again closely resembles a half sine wave pulse. The 

outgoing wave which has had. a longer time to evolve has evolved into 

what can be identified as 5 solitons. The evolution into these 

solitons is shown in the remaining traces. The increased velocity 

with respect to an ion acoustic wave is evident from the trace 

at r = 0.5 cm. The average velocity of the largest ingoing soliton 

is found to be approximately 1 . 17 c* 

s 
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The probe was left in one position and the width of the 
applied pulse varied to determine if the number of solitons was 
dependent upon the app3J.ed pulse width. Fig, 13 shows the received 
signal at .5 cm and at r = 6 cm as the width is varied. For very 
narrow applied pulses a single soliton is observed. As the width 
is increased the first soliton increases in amplitude, speeds up 
and becomes narrower, a second soliton is seen to form. As the 
width is further increased the solitons continue to increase in ampli- 
tude and velocity. The second soliton merges with the first soliton 
and for very wide applied pulses no individual soliton is observable 
though the applied pulse has noticeably steepened. From the data 
shown for r = 0.5 cm an estimate of the product of the ■width of the 
soliton multiplied by the sq.uare root of the maximum amplitude shows 
it to be constant within 10^ for the first 4 traces . 

To determine the velocity amplitude relationship of the 
solitons, use was made of the fact that for very narrow applied 
pulses where a single soliton was found the soliton ait^litude depends 
on the width of the applied pulse. The soliton amplitude can be 
varied from ~ 1^ to 10^ while maintaining essentially a single soli- 
ton. By following the propagation of the pulse the velocity can 
be determined. The propagation of the point of maximima amplitude 
was followed at 1 cm intervals for several soliton amplitudes. These 
results were least squares fit to give the results shown in Fig. l4. 
From this the velocity amplitude relationship for cylindrical soli- 


tons was found to be: 
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V = 


[1 + (1.05 ± .2) 


^1 ]c 
n 'max s 


The measurement of the velocity-amplitude relationship was 
simplified by the fact that over most of the ingoing trajectory the 
con 5 )etition between geometrical growth and damping balanced (see 
Fig. 12). The soliton amplitude was essentially constant over the 
ingoing portion of the trajectory. For a purely ingoing linear 
wave we would expect a (r^/r)^^^ growth to occur. .After accounting 
for this geometrical increase the ion acoustic pulse is calculated 
to have damped by a factor of 5 in propagating 9 cm. Making the 
same correction for the soliton in Fig. 12 we calculate a damping 
of a factor of 1.5 in traveling 9 cm. It must be remembered that 
the (r^/r)^^^ growth should hold only for the linear wave. The 
observed differences may be due to either differences in damping 
or differences from non-geometrical growth. 


C. Discussion 

We can summarize the properties of two dimensional solitons 
be stating that compressive density perturbations evolve into soil- 
tons. The number of solitons is determined by the width and ampli- 
tude of the initial density perturbation. Raref active density 
perturbations do not evolve into solitons. The solitons retain their 
identity after converging (colliding) at the center. All of these 



agree with the well known properties of one dimensional solitons. 

In addition, we find that the soliton width multiplied by the 
square root of the soliton amplitude is approximately constant 
even though the amplitude, width and velocity are functions of 
time. The velocity is greater than the velocity of a corresponding 
one dimensional soliton. Both of these additional properties held 
for three dimensional solitons and the first also holds for one 
dimensional solitons . 

It has been seen from these data that the properties of 
cylindrically symmetric solitons are remarkably similar to these 
of one dimensional solitons . Since no analytic theory exists for 
solutions of the modified KdV equation it was of great interest to 
investigate very carefully the detailed differences and similarities 
that occur . 

It was very difficult to state when the propagating pulse had 
reached a final assymptotic state from which to evaluate properties. 
To overcome this the scaling properties of the ordinary KdV and 
then the modified KdV equation were investigated. This was first con 
sidered by Berezin and Karpman [1967]. 

The ordinary KdV equation is written as: 
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We wish to consider a class of initial conditions of the form 


n(x, o) = “ f (|) 

where ^ is the amplitude of the initial perturbation and L is the 
width. f(A.) is the functional form of the initial perturbation, 
e.g., square pulse, half sine wave pulse, etc. We now transform 
to a new set of variables where the transformation depends upon 
the initial condition according to 


X = 


X 

L 


~ n 
^ ~ 6n/n 



Upon substituting into the EdV equation and simplifying we obtain 


dn 

ht 


+ 


6n A ~ 
— yn 
n L 



1 _A_ 



0 


We have the freedom to choose A such that the coefficient of the 
second term, vanishes-. Upon simplifying we obtain 


Sn Sn 
— + n — + 


St 


^x 




^ , _2 2 en 

= 0 where cr = l — 

n 
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This was obtained using the following transformation 


X = 


X 

L 


~ n 
^ “ 6ri/n 


7r 6n t 
t = — - 
n L 


Under the given transformation all initial conditions are 
transformed to a unit depth -unit width initial condition. Conditions 
with the same numerical values for a evolve identically in the scaled 
coordinate system even though the evolution in the laboratory may be 
quite different. The number of solitons which evolve is determined 
by the numerical value of a. 

There exists a very wide variety of initial conditions for 
which the evolution in this sealed coordinate system is identical. 

The number of solitons in each case must be identical. It comes, 
then, as no real surprise that the dimensionless scaling of the time 
independent Schrodinger equation is the same as the scaling for the 
ordinary KdV equation. The verification of this scaling for planar 
excitations and measurements of the scaling for the cylindrical case 
is of interest and has not been considered before. In order to 
carry this out the deformable screen DP device was used. 

The importance of the dimensionless scaling method is that 
it applies for all time and not only in an as^ymptotic region. The 
known analytic theory provides a connection between the initial con- 
dition and the final assymptotic state where the solitons are 
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spatially separated. The use of the dimensionless scaling method 
allows a direct comparison of soliton formation during the evolu- 
tion from the initial condition. It is only needed to he able to 
identify a point in the evolution process that is easily identi- 
fiable for all initial conditions . This point in the evolution 
must also be able to be measured quite accvirately. The point of 
emergence of the first soliton from the remaining envelope was 
chosen . As can be seen from Fig. 15 the formation of the local 
minimum between the first soliton and the remaining envelope of the 
initial condition occurs quite quickly. It is easily and accurately 
measured and appears for all initial conditions that evolve 
into solitons . 

A large number of data points have been measured for a wide 
variety of initial conditions. The dimensionless time^ , has been 
calculated as has the value of the dimensionless scaling parameter cr. 
These data are plotted in Fig. 16. The points represent both planar 
and cylindrical initial conditions. Both initial conditions are 
uniformly distributed, so no distinction is made between 
them. The points are distinguished by varying values of 

6n 

— I of the initial condition. A systematic variation with — I is 
seen. Larger amplitude initial conditions require a longer time to 
evolve in the scaled coordinates than an initial condition with identi 
cal value of a but a smaller amplitude initial condition. While we 
do not know the analytic form to be expected we do expect the points 
to fall on a smooth curve if the scaling holds. 
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There are several sets of data that were taken for fixed 
initial condition amplitude and varying the width. These are 
shown in Fig, 17 where a differentiation is made between planar 
and cylindrical initial conditions. The same functional dependence • 
is observed for all sets of planar data. The slope in each case is 
nearly identical but there is a consistent deviation with the initial 
value of — The cylindrical solitons^ however, do not display the 
same functional dependence as is evident from the figure. 

Differences do exist for the scaling of the modified KdV 

n 

equation for a cylindrical geometry. The presence of ^ term does 

not affect the transformation used as it will scale the same as the 

— term. 
dT| 

The presence of the ^ term, however, restricts the choice of 
initial conditions , For the ordinary KdV equation a transformation 
of the type 


X = X + a 


t ' = t + a 


leaves the form of the equation invariant. We can then define n(x, o) 
to be an initial condition centered at x = 0 at time t s= 0. Making 
the substitution x = §,t =7] and using the above translation the 
form of the modified KdV equation becomes 
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Sn 7^ 1 n 

w 2Tt’ TTy - ° 


which is no longer the initial equation. 

If we consider the class of initial conditions of the form 


n(x. 


t ) = 
o 


n 


fil) 


the scaling then converts the modified KdV equation into 


3 

5n , Sn . 1 ^ n , n 

at ” ax „ 2 , 5 2t 

2a Sx 


0 


The initial condition becomes a unit width-unit depth initial con- 
dition but since r scales as r/L it is situated at r /l. Figure 

o 

iB shows graphically the effect of the transformation for the ordi- 
nary KdV equation and the modified cylindrical KdV equation. In the 
scaled coordinate system the initial radius of the perturbation 
depends upon the width of the initial condition since the initial 
radius r^, in Debeye lengths, is approximately constant for all data. 
The importance of the ^ term depends upon the initial radius in the 
scaled coordinate system. For this reason we would not expect 
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identical values of ct to exhibit the functional dependence as the 
width is varied to change ^ . This is in agreement with the data 
in Fig. 17 . 

There exists a wide variety of equations that support soliton 

solutions, included in these are attempts to modify the ordinary 

KdV equation to include more realistic plasma effects. Tappert 

[ 1972 ] has included finite ion temperatures. Ott has included 

damping due to ion-neutral collisions [1971] Ott and Sudan 

have treated Landau damping [1970]. 

S charnel [19731 derived a modified KdV equation that allows 

for trapped electrons 


‘ + -h ^ ^ 

where b is a parameter that is the measure of electron trapping. 
Under the transformation 





n 


6n 

n 


we obtain 



~l/2 hn 
- n ' — + 

1 On 

^t 

Sx 

^ ~ J 

2o‘ Bx 
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where 

rSn^'^,2,1/2 
C = [— L ] ' 

Recomputing t and a and assuming b = constant does not significantly 
improve the results. The value for b is determined by the electron 
trapping. Only a weak dependence on — is allowed so the assuit^tion 
of constant b is reasonable. This modified KdV appears insufficient 
to explain the deviations in Fig. I7. 

Tran and Hirt [197^1 have investigated the effect of a two 
species plasma deriving a modified KdV eq.uatiQn 


2 1^ + rn II 

ail as 



The value of 7 is a very sensitive function of the relative con- 
centration and mass ratio of the second species which can arise 
from impTirities in the plasma. A 2^ concentration of hydrogen in 
an argon plasma changes 7 from 1 to .067. The dimensionless time 
and scaling parameter appropriate to the equation are: 

^ ^ ^ a « [7 6n/n 


T 
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Both T and a are functions of 7 . One could, attempt to explain the 
systematic deviation with ^ in Fig. 17 from this. It would not 
be consistent with the manner in which the' experiment was performed 
•however. Much of the data was taken with identical plasma and 
vacuum conditions. Increasing the amplitude and decreasing the 
width to obtain the same value of c- yielded a different dimension- 
less time for the evolution. Several measurements were taken within 
a short time interval. 

If we choose the time sca3j.ng of the ordinary KdV equation 
and search for a functional form of o that gives a smooth curve we 
obtain the results shoira. in Fig . 19- The value for cr is found to be 

a - 

' n 

Wo equation has been found that will yield this scaling. It 
may be necessary to change the scaling of both a and t. 
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IV, CYLmDRICAL ION-ION BEAM INSTABILITY 

A, Experimental Measiirements 
It has already been seen that by applying a DC potential 
between target and driver plasma a steady state radially ingoing 
ion beam can be established (Fig, 4), For a limited range of beam 
velocity a coherent ion-ion beam instability is produced [Hershko- 
witz et al . 197^ ]• This is the result of cylindrical standing waves 
resonating with ingoing and outgoing beams. This experiment differs 
from those in a conventional DP device in that the boundaries play 
a very important role. The inner grounded screen requires vanishing 
potential on its surface. 

The instability is detected by a positively biased Langmuir 

probe which collects electron saturation current. The signal is 

* 

then analyzed using a spectrum analyzer. The power spectrum for 
various beam velocities is shown in Fig, 20. For beam velocities 
greater than those shown the instability abruptly ceases as it does 
for beam velocities slightly less than those shown. The beam velo- 
cities are determined by digitizing the energy analyzer for each 
beam velocity. The digitized traces for a given beam velocity 
are added 64 times to reduce the level of fluctuations due to the 
instability. The averaged traces are then least squares fit to the 


Tektronix 1L5 plug in unit. 
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svan of two Gaussians with the width restricted to he the same. 

A typical fit is shovm in Fig. 21, In. this manner the team 
velocity can he accurately determined even in the region where 
the heam is not resolved. 

The instability manifests itself as a series of discrete 
frequencies that are approximately harmonically related as was 
seen from the power spectrum. A simultaneous comparison of detected 
signals at different positions \ising the two probes showed the 
oscillations to be in phase everywhere. This suggests cylindrical 
standing waves. To verify this the power spectrum was observed 
at several radial positions. The results are shown in Fig. 22, 

The lowest mode amplitude decreases monotonically as the probe is 
moved from the center to the outside. The second mode, however, 
shows a node at r = 7 cm and an investigation of the third mode shows 
2 nodes between the center and the outside. This further enforces 
the identification, of the instability as cylindrical standing waves. 

From Fig. 20 it is apparent that instability frequency depends 
upon beam velocity. A plot- of this is shovm in Fig. 25 where we 
see that the frequency is linearly proportional to heam velocity 
for the lowest 5 modes seen. For the lowest mode we obtain from 
Fig. 25 


(0 5 = .66 v ^ 
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Before attempting to compare these results with theory we must 
determine the distribution function appropriate for the beam. This 
'is accomplished using two energy analyzers; one oriented in the 
radial direction and the other in the azimuthal direction. A beam 
is observable by the radial energy analyzer for all radial positions. 
The azimuthal analyzer can detect a beam in a region near the center 
of the chamber. The size of the region was seen to depend upon the 
energy of the beam. Figure 24 shows a comparison of radial and 
azimuthal energy analyzer traces for a beam of energy 6.0 eV. The 
change in amplitude of the background plasma as the radial analyzer 
approaches the cylindrical screen is due to depletion of the plasma 
in the region between the analyzer and the screen. This depletion 
arises from the blocking of ion trajectories by the energy analyzer. 
Ions that would normally pass through the region occupied by the 
analyzer and reflect from the region near the screen can no longer 
do so. The radial energy analyzer detects a beam for all radii 
while the azimuthal analyzer detects a well identified beam for 
radii less than 4 cm. If we decrease the beam energy the radial 
analyzer shows similar traces. For a beam of 5.0 «V the azimuthal 
analyzer detects a well established beam for radii less than 5 cm 
as can be seen from Fig. 25, For still smaller beam energies which 
appear as barely resolved the azimuthal and radial energy analyzers 
detect similar distribution functions for r k y.O cm. 



When the 'beani is accelerated the "width in energy remains 
constant, Tliis implies that the temp e r at "ur e in velocity space is 
drastically reduced. This allows one to approKxmate the beam by 
a delta function. For the observed distribution function we can 
apprcocimate it by 


f (v) = 6 (v - V ) 

^ 2nv^ ^ 


, 2 2 , 1/2 
V = (v, + ) ' 


where we are in cylindrical coordinate velocity space. 

The beam is represented by a ring of particles of energy 

in velocity space. The particles enter the target plasma with 

energy but the velocities are distributed over an angle. The 

envelope traced out by the angular spread over all points on the 

cylinder represents the boundary over which the ring distribution 

should hold. This, is shown graphically in Fig. 26. 

There are two possible causes of the angular spread of the 

beam as it enters the plasma. The incoming ions have a finite ion 

tenperature of ~ .2 eV. The ratio v /c. represents the tange nt 

Jj X iET* 

of the angle where 20 is the angular spread, where ~ j “~* 


For a two volt beam 



For the system used "we would expect the region to have a radius of 
5 cm. A further increase in the angular spread can arise from 
the passage through and hetween the two fine mesh cylindrical 
screens. If a particle passes through the sheath surrounding one 
of the screen wires it will be attracted or repelled depending 
upon the potential of the wire. Conservation of angular momentum 
a3J-ows only those ions with small impact parameter or velocity to be 
captured by the wire. The particle will follow this perturbed 
orbit in the region between the screens and undergo another 
randomizing passage through the inner cylindrical screen. In both 
of these mechanisms the lower the beam velocity the greater the- 
angular spread will be. 


B. Discussion 

In order to predict the instability and its dependence upon 
beam velocity the Vlasov equation in cylindrical coordinates was used: 


dt 


+ V 


Bf 

hr 




- E 
m r 


.11 


af 

^V„ 



V)251 hf 
^-0 ■ 


We now look at the restrictions on the distribution functions for a 

stationary homogenous equilibriiim by setting =0, e = 0 

Sx 

to verify the restrictions will be compatible with the ring 


distribution function. We obtain 
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■V j 

A 


X — ^ 

0 Sv 


- V 




= 0 




2 2 

The solution to this requires that he a function of v =Jv^ + 


only, i.e. 


fo(^) = f„(v) = ^ 

This restriction holds trivially for the background ion and electron 
Maxwellians . The ring distribution also satisfies this condition 

I 2 2 

since it depends only on § (v - v^) where v * Under 

these conditions we have a stationary homogenous state, with no 
zero order electric field, about which we can perturb. 

We have observed the ring distribution function to exist over 
a region the size of which is determined by the beam energy. The 
instabilities observed typically occur for beam energies of less 
than 2 eV. For these energies the ring distribution function is 
observed to hold in a region that extends to within 5 cm of the 
boundary. 

We can project the ring distribution function onto one axis 
in velocity space to see if it corresponds to an unstable distri- 
bution unlike the shell distribution function which projects to a 



stable distribution function. The ring distribution In cartesian 
coordinates can be represented as: 




1 


2n^B 


6 (v - Vg) 



Integrating out the v dependence we obtain projection onto the 

y 

V axis as 

X 


= i: 




dv 

y 


evaluating we obtain 

for Iv I ■< v_ 

* x' B 

= 0 |vJ>Vg 

This is plotted in Fig. 27. For very large beam velocities FB'(oo/k) 
for low frequency oscillation so the system will be linearly stable. 
For very small v^ we again expect linear stability. For a range 
of Vg we expect the system to be linearly \anstable. 
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To investigate this we perturb the plasma by writing 
f(r> v^ t) = n^f^(v) + Yf 


Then 



+ V . 




V7e can solve this by integrating along the characteristics from 
t = 0 to time t . 


. ox 

V, t) = I -37 


where the integral will follow the characteristics. Initially 
we will work in cartesian coordinates using plane waves as a 
complete set in which to expand. Here the characteristics are 
given by 

X=X-v(t-T) 

OX 

= y - V* ■ 



if we expand §(r, t) in plane waves and look at one mode. 


f(r, t) 




i(ax + by) - iujt 

0 


then substituting we obtain 


^ r* -* , \ ie a,/ j. \ fdi « -i(av’ + bv - (u)s 

X, t) = — #(r, t) ds e ' X y 


1 

X (av + bv ) — r — where s = t - t 

^ X y V $v 


letting t 0 = we have a contribution if > 0 and then 


t^(r, t) = - ^ §(r, t) 


av + bv T ^f 

X y 1 o 

-i(av + bv - w ) V ^v 
^ X y 


To make the connection to cylindrical geometry, we compute 
the first order density fluctuations. 


n(r, t) = J f^ (r, v, t)dv 


substituting and using the identities 
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o 


X = r cos 9 


y =5 r sin ft 


V = V cos 
X 


V = V sin 

y ^ 


a = k cos a 


b = k sin a 


we obtain 


n(r, t) = ~ $(r, t) J vdv J d0 


kv cos^ (0 - g) 1 
kv cos ( p -a ) - tt) V 


defining 





d^. 


cos 0 1 

cos ^ 


3f 
o 


we have 


n(r^ t) = I $(r, t)k^ G(|, f^) 


A specific form was previously chosen for $(r, t). We now choose 
§ to depend on ex according to: 


$ 
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We can recover cylinder f^Inctions for the perturbed density hy using 
the identity 


z (kr) = 1 da, «)) lv((« - tt/ 2) ivJfS _ 


The explicit cylinder function obtained depends upon the path chosen 
to evaluate the integral. For the boundary conditions of the experi- 
ment $ = 0 on the inner cylinder we obtain v = o and lor = where 

Ot is the uth zero of the Bessel f-unction J (p). To obtain the 
oM- o ^ 

dispersion relation we substitute the density 


n(r^, t) = — f J (fcr)k^ G f ) 
^ ^ m o o ^ o 


into Poisson's equation 


V $ = - UiT S e .n . 


which yields 
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D(k, u)) = 1 - E u) G Cr> f 

. po a k 03 ' 
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= 1 “ E 0) r” vdv r^dA . 

. pj *^o kv cos 9 “ u) V ov 


kv cos 0 1 


The terms for background electrons and ions can be reduced 
to derivaties of Fried-Conte functions 
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an explicit form, for the beam contribution due to a delta function 
region velocity space can be calculated for 
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For |Vgj > |~[ we must make an analytic continuation of the 

dispersion relation for real The contribution would be zero 

otherwise. This is analogous to the deformation of the contour used 

to evaluate the Fried-Conte functions for 7 < 0, 

Using the UCSB (University of California at Sata Barbara) 

on-line computer system the dispersion relation with a fluid 

equation contribution for the ion and electron term has been solved 

by iteration. The value of k is held fixed and the iteration is 

performed for a wide range of beam, velocities. This is repeated 

for several beam densities. Figure 28 shows a plot of m vs v and 

b 

7 vs v.„ for the lowest mode with n_ = .3. A maximum in the growth 
B B 

rate is observed for a beam velocity of 1.4 x 10 cm/sec. The ion 
acoustic velocity is 1.5 x lO"^ cm/sec for this case. In the region of 
maximum growth rate the next part of the frequency is proportional to 
Vgk with 



0) = ,61 V k 
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If we express the dependence of the instability frequency 
on beam velocity as 


CD = QV k 

13 


we find the parameter a remains approximately constant as the beam 
density is changed. This is shown in Fig. 50. The experimental 
value of .66 falls well within the observed range of cc. 



V. CONCLUSION 


In a cylindrical geometry compressive ingoing ion acoustic 
pu ls es have been seen to evolve into "cylindrical solitons". The 
propeirbies of these cylindrical solitons have been found to be 
consistent with those of one dimensional solitons (previously 
predicted and observed) and three dimensional solitons (previously 
seen in numerical solutions). The number of solitons is determined 
in a manner similar to the one dimensional case, A detailed investi- 
gation of the evolution of both planar and cylindrical solitons 
showed that differences exist during the evolutionary phase. For 
planar solitons the scaling of the Korteweg-de Vries equation was 
found to be insufficient to explain the observed scaling. However^ 
the scaling has the 'same functional form for all cases. Detai3.ed 
scaling arguments were found to be difficult to apply to the evolu- 
tion of cylindrical solitons . 

An ion-ion beam instability was observed in the presence of 
steady state beams for a limited range of beam velocities . The 
beam was found to be a ring in velocity space. As the beam velocity 
was increased from zero there was an abrupt onset of the instability. 
The instability manifested itself as cylindrical standing waves . 

The frequency of the instability was seen to depend upon 
the beam velocity and an upper 3.imit was found for the beam 



ks 


velocity alDove which the instability was not observed* Numerical 
solution of the dispersion relation appropriate for the standing 
waves yielded good agreement with the observed instability. 
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Figure 1 


Cylindrical double plasma device with two insulated 

concentric plasmas. The inner plasma is produced 

in the right half of the- device-. DC beams- can be 

controlled by varying $ . If $ is pulse a wave 

X3 i3 


can be launched. 
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Figixre 2 


Schematic diagram of the ion energy analyser. Two 
grids and a collector are used. The first grid is 
allowed to float negative. The second grid is 
ramped positive while the collector is "biased at 
-65 volts. Ground shields are placed "between the 
grids snd the collector to prevent capacitive pickup 
between various elements. 
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Figure 2 
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Figure 5 


Current collected by the energy analyzer is displayed 
versus the voltage applied to the retarding grid for 
several different beam energies. The vertical scale 
is arbitrary and the horizcntal scale is 2 volts per 


division. 




Figure 5 
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Figure 4 


Power spectrum of the electron saturation current 
(top) and a differentiated energy analyzer trace 
(bottom). For the bottom trace the vertical scale 
represents the ion distribution function f(v) 
plotted versus the energy to the retarding grid 
(l volt per division). 





Figure 4 
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Figure 5 


A comparison of the ion distribution function as seen 
by an energy analyzer when plotted versus the voltage 
to the retarding grid (top) and when plotted versus 
the ion velocity v (bottom). When plotted versus E 
the density and "apparent" temperature are equal. 

When replotted versus v we see that in fact the beam,, 
density is much less and the beam has effectively 
cooled by its acceleration. 
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Figure 6 


Propagation of a low amplitude ion acoustic wave 
in a fully cylindrical geometry. The vertical 
scale is pertirrhed electron number density 
(arbitrary units) as a function of time (the 
horizontal). The width is 200 M-sec. The top 
trace is taken at r = 9 cm end suceeding traces 


are at 1 cm intervals . 




Figure 6 



Figure 7 


Scheznatic of the modified double plasma device* !Die 


separation screen shown here as a half cylinder is 
continuously deformable (from the exterior) to a plane 
separation screen. The discharge is maintained by 
filaments in both chambers as in a conventional -DP 


device . 
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Figure 8 


Interferometer trace for a planar separation screen. 

The wave amplitude with respect to a fixed phase of 
a reference signal is plotted as a function of distance. 



6h 
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Figure 8 
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Figure 9 


Interferometer trace for a half cylinder separation 
screen with radius 10 cm. • 





Figure 10 


Computed interferometer trace for half cylindrical 
•wave. The values of k , k from the planar case of 
Fig. 8 were used. The amplitude of the appropriate 
Hankel function was computed and plotted versus distance. 
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Figxire 11 Electron density as a function of time at r = 0,5 cm, 
following the application of a half sine wave pulse. 
Signals are labeled "by applied pulse amplitudes. 
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Figure 11 
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Figure 12 Perturbed electron number density as a function of time 
at several radial positions . The dotted traces repre- 
sent linear < 1^) ion acoustic pulses (with ampli- 
tude adjust for comparison). The solid traces represent 
nonlinear pulses propagating, steepening and breaking 
into solitons. 
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Figixre 15 The lefl portion shows the perturbed electron ni^mber 

density detected at r == 0.5 cm versus time. The right 
hand side shows the receii'ed signal versus time at 
r = 6 cm. The traces are labeled by applied pulse 


widths . 



20% div 


74 


0-G74-26 



TIME ( 20 /xsec/div) 


Figure 15 


8% div 




Figure l4 Velocity of a single soliton as a function 
maximum soliton amplitude . 


f the 
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Figure l4 



Figure 15 Perturbed electron number density as a function of 
time at 1 cm intervals. A local minimum is seen to 
form bet-ween- the first soliton and the remaining packet 


8 cm from the screen. 
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Figure 15 
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Figure l6 


6n 


Vli / 

The dimensionless time t = — t/L plotted versus • 


n 

the dimensionless scaling parameter a = 

Points for both planar and half cylindrical Initial 
conditions are plotted and both are distributed uni- 


formly over the graph. 
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Figure l6 



Figure 17 


The dimensionless time is plotted versus the dimen- 
sionless scaling parameter for several sets of data* 
In each set the initial amplitude is held constant 
and the width of the applied pulse is changed to 
vary the value of cr. 
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Figure l8 The effect of the scaling transformation is shown 

graphically for both the ordinary KdV equation and the 
cylindrical KdV equation. For the ordinary KdV 
equation all initial conditions raay be transformed to 
a unit width-unit depth initial condition at the origin. 
The cylindrical KdV equation, however, transforms to 
a unit width-unit depth initial condition located 
at a radius determined by the initial radius in the 
lab and the width of the applied pulse. 
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Figure 19 The dimensionless time is plotted verstis a fit scaling 
parameter. The function form of the scaling parameter 
was chosen to yield a smooth cijrve . 
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Figure 20 Power spectrum of detected electron saturation current 
during instability, showing the first 5 harmonics as 
a function of beam velocity. 
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Figure 21 A comparison of the detected energy analyzer trace 

and a fit to the sum of two Gaussians. This is used 
to determine the beam energy very precisely even 
in the region where beam and plasma are not well 
separated. 
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Figure 22 Pcswer spectrum for the lowest 2 modes shown at 
varying radial positions. 












Frequency versus beam velocity for the lowest 5 

harmonics. No unstable modes are observed for 

5 5 / 

v_ < 0.9 X 10"^ cm/sec or v > 2.2 \ 10 cm/sec. 
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Figure 25 





Figure 24(a) Ion distribution function versus voltage applied t 
the retarding grid fox” the radial energy analyzer. 
Traces are labeled by radial positions. The beam 
is clearly visible foi” all radial positions. 
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Figure 24(b) Ion distribution function versus voltage applied to 
the retarding grid for the azimuthal energy analyzer, 

■ ' Traces are labeled by radial position and beam, energy 
is identical with that of Fig, 24(a), The beam appears 
as well resolved at r = 4 cm. 
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Figure 24(b) 





Figure 25 


Azimuthal energy analyzer traces for a 2.5 eV beam. 
The beam appears • •well resolved at r = 5 cm* 
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Figure 26 Graphical representation of how a radial beam with 
an angular spread of 29 gives rise to a uniform 
' ' density and ring distribution function in a region 

centered about the origin. 
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Figure 27 The projection of the ring distribution function onto 
one axis. An integrable sigularity occurs at v = ± v . 

D 
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Figure 28 Frequency (m) versus beam velocity and growth rate 7 

versus beam velocity for — = .5. In the region of 

i 

maximum growth rate od — ov k where <7 = .6U5. This 
plot arises from solution of the dispersion relation. 
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Figure 29 Plot of o, versus beam density where (X is determined 

from (i) = ov^k in the region of maximum growth rate. 

a 

The values of a are determined from numerical solution 
of the dispersion relation for various values of beam 
density. 
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Figure 50 


Memory allocation of data taking program* 
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Figure 51 


Schematic of location of program control switches . 
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Chip placement diagram for fast buffer memory system. 


Figure 52 
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Figure 53 


Circuit diagram for fast duffer memory system. 
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APPENDIX 

DIGITAL DATA TAKING SYSTEM 
A, General Description 

A versatile computer operated data taking system was con- 
structed using equipment available at the University of Iowa. The 
system consists of a Control Data Corporation (CDC) l6o computer 
and associated peripheral equipment and an interface to the experi 
mental appa,ratus . A Biomat ion 6lOB transient recorder which 
serves as a fast analog to digital converter (ADC), digitizes the 
data and transfers it through the interface to the computer for 
processing. The Biomation 6lOB can digitize 2^6 data points with 
a resolution of 1 part in 64 and store them in its own internal 
memory for later readout by the computer. Data can be digitized 
at rates up to 100 nanoseconds per point. 

The computer equipment presently on line with the system is: 

1) CDC l60 Computer (4096 words of 12 bit memory) 

2) CDC 165 Mag. Tape System (2 tape drives) 

3) CDC 166 Line Printer (250 lines per minute) 

4) CDC 161 Typewriter Station 

In addition, the nuclear physics laboratory has a wide range of 
compatible CDC equipment, including a CDC l6o4 computer system 
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vith Fortran capability. Tapes can be generated by the CDC l6o to 
be processed on the l6o4 system. The CDC l60 computer can be 
replaced with no changes to the interface by a CDC 16O-A computer 
which has options for up to 32 k of memory. 

Interaction of the experimenter with the computer and vice 
versa is accomplished by digital to analog converter (DAC) and by 
means of two rows of switches connected through the interface. 
Programs are generally written to scan the state of the switches. 
Changing the state of a switch will cause the program to branch to 
a predetermined routine , In this manner quite general purpose 
programs can be written. The DAC's are used to drive scope displays 
for the processed data to be displayed upon. 

With the system described here, many useful functions are 
possible. Several of them are: .For low amplitude signals or 

signals in which the noise is comparable to or greater than the 
signal, by digital addition of many passes the relative noise level 
is reduced and a clean signal is recovered. Even in the cases where 
the noise level is quite small, repeated addition of successive 
passes gives better resolution of the wave form. Precise measure- 
ments in time of flight experiments is possible with digital data. 
After repeated additions, the maximum (or minimimi) is lisually 
identifiable as one channel. For data for which an expected form 
is known a computer fit to the data is possible. An exaniple of 
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“this is seen in Fig* 21. Through "the use of digital to analog 
converters, voltage levels from the computer can he fed to 
experimental apparatus, e.g. voltage controlled oscillators, 
remote programmable power supplies , and wave form generators , etc . 
With two Biomat ions correlation function measurements become possible 
With the existing equipment these and many more measurements 
are possible, limited only by the time resolution of the analog to 
digital converter (Biomat ion 610B), 

The initial work on the interface from the Biomation to the 
CDC 160 computer was done by Dr. Stephen Wender [1975]* The 
wiring to convert the CDC logic levels (-18,0 volts) to the standard 
TTL logic levels (0,5 volts) and the logic circuitry necessary to 
decode the computer selection of devices was done first. Two rows 
of 12 switches each (used for program control) and circuitry 
for inputing from the Biomation were also completed by Dr, Wender. 

In his initial design there were 8 input ports that could be used. 

Two were used for program control switches and one was used for 
the initial Biomation 6IOB. A second Biomation has recently been 
added to the system. The four remaining parts can be used for any 
device with digital output (e.g., a digital voltmeter). 

The author has extended the previous work to include provisions 
for 8 output ports. One digital to analog converter (DAC) can be 
connected to each port. This allowed for computer feedback of 
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of voltage signals to the experiment. The most standard use for 
the DAC is for scope displays of the digital data. They are also 
used to drive an x-y plotter for graphical copies of the processed 
data , 

A fast buffer memory system was also designed and constructed 
to increase the transfer rate between the Biomation and the CDC l60 
con 5 )uter. The Biomation can output the digital data at a rate of 
either 2 microseconds per xTOrd or at a rate greater than 512 micro- 
seconds per word. The CDC l6o computer can input information at a 
rate of 4.8 microseconds per word. To obtain maximum transfer rate 
the Biomation dumps its memory into a buffer memory at 2 microseconds 
per word and when it is done the computer reads the buffer memory 
at 4.8 microseconds per word. This decreases the time for transfer 
of one data pass from a minimimn of .13 seconds to .0015 seconds. 

B. Program Description 

1 . Program 

Data-Tal.e, D-A, with Autoscale 

2. Purpose 

General purpose data taking program. Signal averaging by 
simple addition up to 64 times. Provisions for storing and viewing 
up to 12 data sets are included; data can then be dumped on magnetic 
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tape with an identifying header of up to 256 symhols. Program 
control is hy means of the switch box. 

5. Memory Usage 

Locations OOOOg to location lOOOg are used for program 
storage. Locations lOOOg to ikoo^ are used for scratch storage. 
Locations lUoOg to 2000g are termed working storage. Locations 
2000g to TTTTq s-re divided into 12 storage areas of 400g = ^^^10 
words (see Fig, 50 ). 

4. Program Control Switches 

Two rows of twelve switches are used for program control . 
(Fig. 31 ). Each row represents one computer word and is grouped 
into four sets of three switches . Each set of switches rei^resents 
one octal digit. Thus, any rovr of switches can be set to represent 
any number from OOOOg to 7777g. The three switches that represent 
an octal digit are a binary representation of that digit (e.g., 
on-off-on represents 5). A switch is given a 1 value by setting 
and leaving it toward the top. 

Switches 1-6 on the bottom row are used for program control. 
The remaining I 8 switches are used to feed information to the pro- 
gram as it executes. Switches J-±2 of the bottom row are used to 
tell how many times a signal should be averaged. They can be set 

from 1 to 77g = 
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Switch 1 on the hottom row is used to increase the versa- 
tility of* the prograni. The role of switches 2-6 depends crucially 
on what position switch 1 is set to- For switch 1 = 0: 

Switch 2: Causes retake of data when set to on position and 

then returned to off. Averaging will take place and a return to 
the display mode will he made. In the display mode the last data 
set taken is displayed. This is stored in the working storage 
area of Fig. 1. 

Switch 3 : When in display mode setting this switch on then 

off again causes the displayed data to be multiplied by a factor of 
Switch 4: When in display mode setting this s'witch on then 

off again causes the displayed data to be divided by a factor of 2. 

Switch This is used to store data into one of 12 avail- 
ablve storage areas . The storage area is determined by the top row 
of switches. Switches 1-^ of the top row are used for this. The 
computer treats the remaining 8 switches as if they are zero, 
regardless of what they are set to be. The top row is then treated 
as the first word address to which the data is to be moved. This 
can take on any value from 2000g to 7^00^. Switches 1-5 are used 
to determine the first octal digit of the location, e.g,, 2000, 
4000, etc. The fourth switch is used to determine if the data is 
stored at xOOOg if the switch is zero, or at xUoOg if the switch is 
set.- Caution. It is possible to destroy the program if a store is 


iT> 

c. • 


made to location OOOOg or 0400gl 
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Switch 6 ; This is losed to view the 12 storage areas . Switches 
1“4 of the top row select the appropriate storage area. They can he 
changed while switch 6 is set. The remaining 8 switches on the top 
row are used to zero one location for display purposes. This allows 
comparisons of different sets of data. These 8 switches can take 
on values from 0 to 577 q corresponding to the 256 channels of the 
Biomat ion. 

There are two 8 hit D-A converters available to the system. 

One of these is connected to the lower 8 hits of the 12 hit computer 

V7ord while the other is connected to the upper 8 hits. Depending 

upon the number of times averaged, it may he more convenient to 

use one of these over the other. Provisions for this have been 

made in the program and a choice is implemented when switch 1 of 

the lower row is set to 1. Switch 1 set to 1: 

Switch 2 ; Used to dump data onto magnetic tape. When 

switch 1 and switch 2 are thrown in that order, the input light on 

the typewriter station should come on. A header of up to 256 

characters can then he typed in. The symbol "/" is not allowed as 

a valid s^rmhol. At the end of the header a is typed and this 

initiates the dumping of the header and the 12 storage areas in the 

computer onto magnetic tape. A return is then made to the display 

mode. Switches 1 and 2 should he returned to the zero position 

* 

before typing in the header. 


Tape to he written on must he loaded on a tape drive with unit 
number set to 1, 
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Switch 3: This selects the D-A converter that is connected 

to the upper 8 hits. The switch near the oscilloscope must also 
he in the correct position. To change from, one to the other set the 
switch near the oscilloscope to the desired D-A converter and then 
flip switch 1 followed hy either 5 or 4. Return both switches to 
off and the display should work. 

Switch 4 : This selects the D-A that is connected to the 

lower 8 hits. It is selected as described in the previous section. 

5. To Run 

Load program, clear and run fr*om P := 0100. To initialize 
Biomation remove trigger signal, set to external trigger, clear 
computer and run from 0100. Light on left side of console should 
display static "IN", Remove from run, then clear and i*uh from 0100. 
This assures that the first word input from the Biomation will he 
recognized hy the computer as the first data point. It may he 
necessary to repeat this later if display indicates improper trigger. 
Connect trigger to Biomation and proceed to take data. If problems 
occur, run from 100 to take data or from 153 to go directly to 
display mode. (Note: If it is run from 0100, the first thing the 
program does is try to input data. ) 

6, Program 


l60 Data Plot 
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7* Purpose 

To plot data that has been stored on -siagnetlc tape, 

8. Memory Useige 

Locations 0000 to l400 are used for program storage \7hile 
locations l400 to 7777 used for storing data read in from 
magnetic tape, 

9* Use 

After loading program into computer run from location 0100^ 
all switches on the svzitch box must be set to zero. The top row 
of switches is used to transfer infomnation to the computer while i 
the bottom row is used to control the program. The switches in 
the bottom row will be labeled from 1-12 staring from the left. 

Switch 1 : The left most switch is used to input data from 

the tape. Flipping the switch momentarily inputs one file from, the 
tape and should type out the associated header. If switch is not 
returned to zero position, the tape will continue to read in. 

Switch 2 : Searches forward one file mark on the data tape. 

Switch 5 : Searches backward one file mark on the data tape. 

Switch 4 : Displays the data on the oscilloscope. It is 
used to display several sequential sets of data. The sets to be 
displayed are specified by the top row of switches . The left four 
bits specify the starting address. The program treats the left four 
bits as the uppermost hits of a word with ail other bits set to 0. 



The numiber of consecutive sets of data to be display is specified 
by the lowest four bits of the top row of switches . 

Switch Plots the data on the x-y plotter. The data to 
be plotted is specified in the same way as for displaying with 
switch 4. The upper four bits determine the first data set to be 
plotted while the lower four bits determine the number of plots 
to be put on the same page. 

Switch 6 ; Plots a square to help in the alignment of the 
x-y plotter. The square plotted has the maximum range of x and 
y values possible and thus serves as a border for the paper. 

Switch 7 : Divides all of memory by 2. This is used in 

scaling to the correct height. 

Switch 8 ; Multiplies all of memory by 2 in case the data i 
too small. 

Switch 9 1 Displays locations l400 to 2000 with the channel 
specified by the lower 8 bits of the top row set to zero. This 
can be used to follow the evolution of a waveform. 

Switch 10 : Used to exchange certain spectrum that are 

specified. The contents of locations l400 to 2000 are exchanged 
with the data starting at the location specified by the left 4 
bits of the top row. 

Switch 11 ; Adds the amount specified by the lower 8 bits 
of the top row to the spectrum stored starting in location l400. 
This is used to adjust the heights of different spectrum. 
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S witch 12 t Displays both the contents of locations l400 
to 2000 hut also the spectrum specified, by the left 4 bits of the 
top ravT^ This is used to compare the present data set to which 
one is adding set amount to the previously adjusted data set. With 
switch 12 set, any other switch can also be used. (Note: With 

switch 12 and switch 4 there is a slight problem with triggering 
of the oscilloscope. It is necessary to expand the trace with the 
time base slightly in order to get all of the traces displayed 
properly. ) 

C. High Speed Buffer Memory 

In order to fully utilize the high speed characteristics of 
the Biomation olOB, it is necessary to build a high speed buffer 
memory. With this memory it is possible for the Biomation to read out 
its zoemory at 2 microseconds per word and then for the l 60 con^zuter 
to read the infoimiation at its highest possible operating speed. 

This fast buffer memory is built on a card that fits into 
the original interface built by Dr. Stephen Wender. It makes use 
of the signals available from the interface. Figure 32 shows a 
diagram of the chip placement on the card. Figure 55 is a schematic 
circuit diagram with a list of the integrated circuits used. 

The sequence of events that occur is described below. The 
computer executes a 7500 5400 instruction. When this happens, 
the 54 oO a line from the computer makes a 0 to 1 transition. This 
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causes the output of IC 9 A to inalte a 1 to 0 transition. This 
fires the one-shot ICI 9 . ‘The Q output of this is xised to arm 
the Biomation. The Q is fed into nand gate IClOA, The output 
of this is fed to the clear pin of the 4 hit counters IC15-14. 

After the Biomation is armed, it -will be triggered and the flag 
will go high. The flag is inverted by IC 9 C and fed to the write 
enable memory chips where it causes the information putoiit "by 
the Biomation to be written into the memory chips. This sig- 
nal is also nanded through ICllA to the one-shot ICI 7 . The firing 
of the one-shot advances the counters, so the next word from the 
Biomation will be read into the next memory location. The flag 
from the Bicmation is also inverted twice by IC 9 D and E and fed to 
the word command on the Biomation. This causes the flag to be 
removed until the next word is available . After executing the 7500 
5400 sequence a 7^00 (input to A) is executed. The input request 
line from the computer makes a 1 to 0 transition. This, combined 
with the 54 oO A line from the interface, is fed through IC12A and B 
and through IC 9 F to one input of ICIOB. The other input of the 
nand gate IClOB comes ffom an 8 input nand gate where the inputs 
are hooked to the 8 address lines of the binary counters. It fires 
when the Biomation has read out all 256 channels. The output of- 
ICIOB is fed to the information ready and the computer resumes 
operation. 
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A similar procedure is followed to read the information from 
the memory chips to the computer. The computer executes a 7500 
5401 which causes the 5^01 A to make a 0 to 1 transition. This is 
fed through IC9B and fires the one-shot IC20. This is then fed to 
the clear inputs of the counters. The input request is inverted 
by IC12A nanded with 5^01 A by IC12C. The output of this is then 
fed back to the information ready. It also fires the one-shot 
ICl6. This is nanded throi;igh ICllA to fire the one-shot ICI 7 
The purpose of ICl6 is to put in a delay before ICI 7 fires to 
advance the counters. Input continues at the maximum rate 
possible for the computer. 


Integrated Circuits Used 


No. 


Description 

Symbol 

6 

74 s 206 

256 x 1 RAM Memory 

ICl -6 

2 

54195 

4 Bit Coimters 

IC13-3.4 

2 

7437 

Buffered Hand 

IC7-8 

5 

74121 

Monostable Multi- 
vibrators 

ICI 6-20 

1 

7450 

8 Input Nand Gate 

ICI 5 

1 

74o4 

Hex Inverter 

IC9 

3 

7400 

Quad 2 Input Hands 

icio-12 


Note: The 5^195 could be a 7^195. 

External resistors of 10 K are hooked between the 
outputs of the memory chips and +5 volts. These are described in 
the literature for the 7 ts 206 . 
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Observations of Ion-Acoustic Cylindrical Solitons 

Noah Hershkowitz and Thomas Romesser* 

Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52242 

(Received 24 January 1974) 

Cylindrical solitons are seen to evolve from compressive cylindrical pulses in a coHi- 
sionleas plasma. The properties of these solitons are found to be consistent with the 
known properties of one- and three-dimensional solitons. 


Soliton solutions are now well known for at 
least seven distinct -dimensional wave sys- 
tems, In particular, the soliton solutions of the 
Korteweg-de Vries (KdV) equation have been ex- 
tensively studied both theoretically and experi- 
mentally during the last decade and KdV is now 
known to approximately describe many systems 
which include nonlinear and dispersive effects.^ 
Washimi and Taniuti^ have shown that slightly 
nonlinear one-dimensional ion acoustic waves in 
collisionless plasmas with cold ions are described 
by KdV. In a recent Letter, Maxon and Viecelli,® 
following the procedure of Ref. 2, have derived a 
modified KdV equation for spherically symmetric 
ingoing waves. In this Letter we present experi- 
mental observations of cylindrical solitons in a 
coUisionless plasma. 

Oie- dimensional solitons have several distin- 
guishing characteristics.**’^^ Among them are 
the following: (1) Arbitrary positive (compres- 
sive) density perturbations evolve after sufficient 
time into a superposition of spatially separated 
solitons (solitary pulses). (2) The number and 
amplitude of the solitons is determined by the so- 
lution of an appropriate time -independent Schrd- 
linger equation with a potential well that is pro- 
aortional to the initial spatial density perturba- 
aon. One soliton is formed for each bound state 
vith soliton amplitude proportional to the energy 
jigenvalues, (3) The soliton velocity is given by 
1 +-j(6«/«)]c^, where tn/n is the maximum den- 
sity perturbation of each soliton and is the ion 
Lcoustic velocity. (4) The spatial widths are pro- 
>ortional to {tn/n) which implies that the 
iroduct of the square root of the maximum soli- 
on amplitude multiplied by the width is a con- 
tant. (5) Solitons retain their identity upon col- 


lision with other solitons. 

All of these properties have recently been veri- 
fied experimentally with collisionless plasmas. 
Linear double-plasma (DP) devices were used by 
Ikezi, Taylor, and Baker*^ to verify all but the 
second property, and by Hershkowitz, Romesser, 
and Montgomery^^ to verify the connection with 
the underlying SchrfSdinger equation [property 
(2)]. Cohn and MacKenzie^ investigated solitons 
resulting from very large density perturbations 
produced by photo ionization, A summary of much 
of the experimental evidence has been given by 
Ikezi. 

In the first woric, which considers solitons of 
dimensionality greater than one, Maxon and Vie- 
celli^ have numerically determined that spherical 
solitons have the following four properties. First, 
an ingoing soliton increases in amplitude while 
decreasing in width, thus retaining its identity 
as a single soliton. Second, the product of the 
square root of the maximum soliton amplitude 
multiplied by the width is a constant. Third, a 
small residue develops and moves inward behind 
the soliton, taking up a measurable percentage 
of the total momentum; and fourth, the soliton 
velocity is somewhat greater than the velocity of 
a corresponding one- dimensional soliton. 

In this Letter we present data showing that cyl- 
indrical solitonlike objects exist and that their 
properties are consistent with those of one- and 
three-dimensional solitons. These results are 
to our knowledge the first experimental evidence 
for solitons of dimensionality greater than 1. 

Experiments were carried out using a cylindri- 
cal DP device which had previously been used to 
study the ion-ion beam instability of cylindrical 
beams and background plasma.^® Two concentric 
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cylindrical plasmas (lei^th 30 cm) are separated 
by two closely spaced, fine-mesh, concentric 
cylindrical screens with inner screen diameter 
equal to 20 cm. The outer screen is negatively 
biased to prevent the flow of electrons between 
the plasmas, and the inner screen is grounded. 
The ion density was approximately 10® cm'® and 
the ion and electron temperatures were approxi- 
mately 0.2 and 3 eV, respectively. Positive half-‘ 
sine-wave pulses are applied to the outer plasma 
to launch cylindrical density perturbations in the 
inner plasma. 

Signals are detected by a positively biased 
Langmuir probe which has variable radial posi- 
tion. No azimuthal dependence was observed. 
Figure 1(a) shows the perturbed electron number 
density as a function of time at several radial 
positions for both large and small initial density 
perturbations. For the small -amplitude pulse at 
r = 9 cm we can identify an irgoir^ pulse, which 
is quite similar to the applied pulse, followed at 
a later time by a similar outgoing pulse that has 
propagated from the opposite side and through 
the center. As the probe is moved closer to the 
center the ii^oing and outgoii^ pulses approach 
each other, mergii^ at the center. For the large- 
amplitude compressive pulse at r = 9 cm the ingo- 



FIG. 1. (a) Perturbed electroa number density as a 
function of time at several radial positions. Upper 
traces, linear (dn/n< 1%) ion acoustic pulses (with am- 
plitude adjusted for comparison). Lower traces, non- 
linear pulses propagatit^, steepening, and breakii^ in- 
to solitons. (The received signals are digitized and 
stored on magnetic tape for later analysis. This is the 
cause of the observed steplike structure.) (b) perturbed 
electron number density detected at r=0.5 cm labeled 
by the applied pulse widths, (c) Received signals at r 
= 6 cm labeled by applied pulse width. 


ing pulse is seen to be similar to the applied 
pulse, but three solitons can be identified in the 
outgoing pulse. The traces at other radial posi- 
tions indicate how the initial density perturbation 
evolves into the solitons. The increased velocity 
of the first two solitons compared to the ion 
acoustic velocity is evident. Once formed, the 
largest soUton is seen to be much narrower than 
the applied pulse. We find that the average veloc- 
ity of the largest ingoing soliton is approximately 
1.1 7 C 3 . The application of a negative (rarefac- 
tive) density perturbation is not found to evolve 
into solitons. 

As in the one- dimensional DP device, the maxi- 
mum applied voltage is limited by the electron 
temperature/® For applied voltages larger than 
the electron temperature (here »3 eV), particle 
bursts (pseudowaves) are detected. Data were 
taken with the largest initial density perturba- 
that could be obtained without launching pseudo- 
waves. The width of the applied pulse was then 
varied to determine how the soliton number de- 
pended on the initial density perturbation (see 
Fig. 1). This procedure was identical to that fol- 
lowed in our earlier measurement (Ref. 13). 
Figure 1(b) shows how the signal received, r 
= 0.5 cm, depends on applied pulse width. Widens 
ing the applied pulse results in increased ampli- 
tude, decreased width, and increased velocity in 
the received signal. In the top trace we see one 
well-defined soUton. In the second trace the 
first soliton has grown and narrowed and a sec- 
ond soliton is apparent. In the third trace the 
first two solitons have grown, narrowed, and 
speeded up. A third soliton is barely apparent. 
We find that the square root of the maximum am- 
plitude multiplied by the width is constant to with- 
in 10% for the first four traces. In the fourth 
trace the third soliton is seen to grow as well. 

For further increases in width the solitons no 
longer have sufficient time to separate from the 
initial perturbation. Figure 1(c) shows the sig- 
nals corresponding to the same six applied pulse 
widths as detected at r = 6 cm. 

For small applied pulse width it is possible to 
launch single solitons whose amplitude depends 
on the applied pulse width. The amplitude of an 
incomir^ soUton was found to be approximately 
constant over much of its trajectory as a result 
of a competition between damping and geometric 
growth (see Fig. 1). This simplified the deter- 
mination of the amplitude dependence of the ve- 
locity. Figure 2 shows the soliton velocities, de- 
termined from individual soliton trajectories, 
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FIG. 2. Velocity of single solitons as a function of 
the maximum eoliton amplitude. 


versus soliton amplitude. The best least-squares 
fit to these data gives u =[l 4a(C«/n)]c^, where a 
= 1.05±0.20. This is faster than a corresponding 
one-dimensional soliton. 

This experiment differs from an idealized one 
in at least two respects. First, as in one-dimen- 
sional experiments, damping is present. After 
accounting for a geometric increase which goes 
like the ion acoustic pulse is seen to 

damp by about a factor of 3 in propagating 9 cm, 
and the soliton damps by a factor of 1.5. In the 
absence of damping, the geometrical growth fac- 
tor {rjr) should only be expected to hold for lin- 
ear ion acoustic pulses. Second, the density per- 
turbation is not found to diverge at the center. 

We attribute this result in part to broadening of 
the received pulse as a result of variations in 
the radius of the cylindrical screens of the order 
of 0.3 cm and to finite size of the probe and of an 
insulating glass cylinder (0,5 cm) which covers 
all but the last centimeter of the probe. 

We summarize the measured properties of two- 
dimensional solitons. Compressive density per- 
turbations evolve into solitons. The number of 
the solitons is determined by the width and ampli- 
tude of the applied pulse. Rarefactive perturba- 
tions do not evolve into solitons. The solitons 
retain their identity after converging (colliding) 
at the center. All of these are well-known prop- 
erties of owe -dimensional solitons. In addition 
we find that the soliton width multiplied by the 
square root of the maximum soliton amplitude is 
approximately constant even though the ampli- 


tude, width, and velocity are functions of time. 
The velocity is somewhat greater than the veloci- 
ty of a corresponding one-dimensional soliton. 
Both of these properties of three-dimensional 
solitons and the first holds for one- dimensional 
solitons as well. 

Maxon has recently derived a modified KdV 
equation for cylindrical solitons,^’ Detailed com- 
parison witii numerical solutions of this equation 
will be presented in a later publication. Attempts 
will be made to compare the amplitude, width, 
and propagation speed as a function of time with 
numerical solutions of a modified KdV equation 
for cylindrical solitons which includes damping 
when such results become available. 

We thank S. Maxon for helpful discussions and 
Alfred Scheller for construction of much of the 
apparatus. 
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I'he ion ion instability is studied in a cylindrical double~plasma device. IjOw frequency 
cylindrical standing waves are found which are one-dimensionai in character with fre- 
quency proportional to beam velocity. An approximate dispersion relation for the cylin- 
drical standing waves is derived. 


Double-plasma (DP) devices have recently been 
used to study the incoherent ion-ion two-beam in- 
stability produced by one-dimensional beams in 
unmagnetized plasmas.^"® This turbulence is 
found to be /^?ree-dimensional in character in 
agreement with the linear theory of the ion-ion 
instability^’"^ which predicts that the instability 
depends on k'v^, where ^ is the pZawe-wave prop- 
agation vector and Vj is the beam velocity. Al- 
though there is an upper limit on beyond which 
the o«c-dimensional ion-ion instability will not 
grow, there always will be oblique directions for 
which the projection of v^ on k will give growing 
modes in three dimensions. Means et wZ.^^have 
recently argued that the observation of turbu- 
lence in experiments with owe -dimensional elec- 
trostatic ion acoustic shocks depends fundamen- 
tally on this Z/zrec-dimensional property of the 
instability. 

In this Letter we report the production of a 
coherent ion-ion instability which is essentially 
one -dimensional in character. This has been ac- 
complished by generating cylindrical standixig 
waves which are resonant with ingoing and out- 
going beams in a cylindrical DP device at the 
University of Iowa. This device differs from con- 
vential DP devices^ in that the cylindrical bound- 
ary of the plasma plays a dominant role. A 
grounded cylindrical screen through which the 
beam in injected, serves as a well-defined bound- 
ary condition (vanishing potential) for the stand- 
ing waves. In conventional devices^ "= the dimen- 
sions were such that wave and particle phenom- 
ena were not significantly affected by the pres- 
ence of boundaries. The dimensions of those de- 
vices were lur^e compared to the ion charge-ex- 
change length, the e folding distance for ion- 
acoustic waves, and all wavelengths of interest. 

In the cylindrical DP device described here the 
diameter of the plasma is comparable to these 
lengths. 

A description of these standing waves is de- 
rived from the Vlasov equation in a cylindrical 


geometry. The dependence of the frequencies <j 
of the instabilities on the beam velocity is shown 
to be remarkably similar, but not identical, to 
the results for a strictly one-dimensional ion- 
ion instability. 

The cylindrical DP device, which has recently 
been used to study cylindrical solitons,^^ is shown 
in Fig, 1. Two concentric argon plasmas are 
separated by an outer negatively biased screen 
and an inner screen held at ground potential. 
Plasma within the inner cylinder is produced in 
an adjacent connected chamber. Typical operat- 
ing parameters were electron temperature 

eV, ion temperature T. ^0.1 to G.2 eV, aver- 
age plasma density ~ 10® to 10® cm'^, and pres- 
sure sj 2x10"'* Torr. A steady-state radially in- 
going cylindiucal beam is formed by raising the 
plasma potential in the outer cylinder. Beam 
density ratios Tj=n^/n^ are controlled by varying 
the concentric discharges. The ion charge-ex- 
change length was greater than the inner-cylin- 
der radius. 

Energy distribution functions in the inner plas- 
ma (region A in Fig. 1) are measured with two 
energy analyzers with depth-to-area ratios great - 



FIG. 1. Cylindrical DP device with two insulated con- 
centric plasmas. The inner plasma is produced in the 
right half of the device. Beam energy is controlled by 
varying 4>£. filament supply voltage; filament- 
to-wall voltage; , applied bias voltage, 2’^ = 1.0 eV; 
T^S.0,2 eV; A'f| = 10^-10^ cm'^; pressure, 2xiO"*Torr. 
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FIG. 2. Typical radial and azimuthal energy-analyzer 
traces at 5 and 8 cm are shown. Similar traces are ob- 
served for radii less than 6 cm indicating the presence 
of a ring distribution function in velocity space. The 
absence of an azimuthal beam at 8 cm is apparent. 
Azimuthal traces are found to be independent of 'P . 
Broadening of these traces is instrumental . 


er than 1, whose radial coordinates can be va- 
ried. One energy analyzer measures the energy 
distribution function of particles with velocities 
in the radial direction, another observes the 
energy in the <p direction. A third measures out- 
going energy distribution functions in the outer 
plasma (region B in Fig. 1). Typical energy- 
analyzer traces are shown in Fig. 2. Near the 
inner screen the beam is seen to be radial. For 
radii less than 7 cm we find approximately equal 
radial and azimuthal beam components, indicat- 
ing that the beam forms a ring in velocity space. 
The spatial region over which the ring distribu- 
tion function exists is determined by the beam 
velocity, and the separation and mesh size of the 
two concentric screens. The region increases 
for smaller energies. The presence of a ring 
distribution function rather than a purely radial 
beam results in a uniform beam density (?-<7 
cm) with no steady-state electric field. This 
facilitates a theoretical description of the in- 
stability. 



( a ) 



BEAM VELOCITY IN UNITS 10^ cm /sec 
(bl 

FIG. 3. (a) Power spectra showing the first three 
harmonics as a function of beam velocity . (b) Frequen- 
cy versus beam velocity for the three lowest harmon- 
ics. No growing modes are observed for <0.9 xio^ 
cm/sec and for Vj,>2.2xio® cm/sec. 


Instabilities are detected by positively biased 
Langmuir probes oriented in the axial direction 
and by the energy analyzers. One Langmuir 
probe is variable in the radial direction and the 
other is variable in the <f> and axial directions. 
Langmuir probes indicate that relatively uni- 
form background plasma and beam densities are 
achieved in this device. Langmuir-probe mea- 
surements within the inner cylinder (region A in 
Fig. 1) showed no or ^ dependence (except near 
the cylinder ends). 

A comparison of signals simultaneously ob- 
tained at different positions showed that the in- 
stability was in phase throughout the inner plas- 
ma, demonstrating that standing waves were pro- 
duced. Therefore it is meaningful to consider 
power spectra. Typical instability power spec- 
tra as a function of beam velocity are shown in 
Fig. 3(a). For high beam velocities we observe 
no histabiiny. The onset of the instability is 
seen for a beam energy of approximately 1.0 eV. 
For beam energies less than approximately 0.2 
eV the instability disappears. As the plasma po- 
tential of the outer plasma is made less than the 
potential of the inside plasma we observe inco- 
herent instabilities between the outer cylindrical 
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screen and the cuter walls of the chamber (re- 
gion B in Fig. 1). 

Measurements of the power spectra in the ra- 
dial direction show that the three lowest frequen- 
cies correspond to resonant modes with no nodes, 
One cylindrical node, and two cylindrical nodes, 
respectively, further demonstrating cylindrical 
standing waves. Figure 3(b) shows that w is 
roughly proportional to Uj, for the first three har- 


monics. In addition, the second and third har- 
monics have approximately 2 and 3 times the fre- 
quency of the first. This suggests that = 
with a constant, describes all three modes" with 
/?„ varying approximately with the mode number. 

These results can be understood by consider- 
ing the linear Vlasov equation in cylindrical geo- 
metry. The Vlasov equation can be written for 
one component (ions or electrons) in the follow- 
ing way: 


3 / df 

'\'V ^ 

9/ ^ dy yd(p 




(1) 


For a stationary homogeneous equilibrium distribution /„(r, u,) with E =0. we find ‘(n. afjav 

- n =0. This means that the zeroth-order distribution function depends only on the magnitade 

e velocity, foU>)~fd{v^ +v^ ) ^ ]; i.e., the distribution is concentric in velocity space and the 
beam must be a ring. Experimentally, we find such a distribution function extending from the center 
to within 3 cm of the inner screen for beam energies less than 2 eV, 

Perturbing the plasma, we write /(?, v, t) = +/,(?, v, t), where /.(f, f, t) is given by 

/i(f, V, V*(r„(r, V, I - f„), <„) • 3/„(n)/8v, (2) 

with * being the electrostatic potential. The time integral is taken along the straight-line orbits of 
the unperturbed state In evaluating the time integral we consider the plane-propagating waves $(r. t) 
-■|.<,exp[i(mr v6y _ u,i)J = *„ exp{t[*r cos(<(> - a) - uif]}. Then the integral in Eq. (2) can be performed to 


/,(r,v,/)=^^(r A... fe^cos(e-or) 3^^ 

m * kv cos(6 - oc) - (o vdv " (3) 

In writing Eq, (3) and the potential we have used the following definitions: 
x = rcos<^, y=r sirup, v^ = vcosO, v^ = vsin6, a =1? cos or, b=ksina. 

When calculating the density from Eq. (3) by integrating over n and 6, we note that the dependence of 

cos(e- a) on q; can be suppressed because e is integrated over all angles. Thus we obtain for the 
density 

With 


® cos0-u)/fei; vbv 


(5) 


In ordei to express the density in cylinder functions we let depend on n, 4 . '<,) = (* /xlexnhvo. 

- w/2), and integrate over or. A change of the variable or to a; = A + <^ and an appropriate extension of 
the limits of the integral to infinity produces the integral representation of the HanLl or Bessel func 
tions Z^{kr)eyc^{iv<p), Thus we can write 


( 6 ) 


«(r, /) = {e/m)^,,Z^^ikr) exp(ii/<p)G . 

The densities have to be inserted into Poisson’s equation, 

We note that Z^(fer)exp(?i'i;fl) is an eigenfunction of the Laplacian. This shows that cylinder functions 
are eigenfunctions of the beam-plasma system. For the boundary conditions of the experiment ^ = 0 
on boundary, we obtain directly u =0, kR = where is the pth zero of the Bessel function J (p) 
It IS important to notice that these are the only nonsingular solutions and they represent standing ° 
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waves. The dispersion relation is found from Eqs. (7), (6), and (5): 

D(k, «) = 1 = 1 


a? V 3r 


= 0 . 


( 8 ) 


Note that = Og ^ a plane~y/3.ve propagation vector. 

We consider now a background plasma with Maxwellian ion and electron distributions with T^» Tj. 

In addition we have a radial influx of ions of velocity which can approximately be described by 
f^{v)-{nf/2Tiv^)5{v -v^). From the Penrose criterion it follows that for small the system is stable 
as well as for very large Vj. Both effects have been observed experimentally [compare Fig. 3(b)]. 

From the assumed distribution functions we find for very small ion beam density (t?« 1) and for Uj 

0 )=+feUj[l + exp(2irz‘ /3 ) c/ ^ 

For the case of resonance (t'j, = C 5 ) we find w =+A’t^^l + exp(2ve'/5)]. The experimental proportional- 

ity of o> with is evident from Fig. 3(b) which corresponds to tj«Q.2. 

In conclusion we have shown that low-frequency cylindrical standing waves which depend only on r 
are produced by the ion-ion beam instability in a cylindrical DP device. We have shown that their fre- 
quency is roughly proportional to the beam velocity. We believe that this is the first time that a coher- 
ent ion- ion beam instability has been observed. We ascribe this to the high symmetry of the experi- 
ment which substantially reduces the off-axis modes. We have derived a dispersion relation for the 
coherent cylindrical- standing- wave instability from the appropriate Vlasov equations for the cylindri- 
cal geometry. This predicts unstable standing waves with proportional to i\ rather-^than travel- 
ing waves. 

We are indebted to David Montgomery for advice, assistance, and enlightening discussions. We are 
grateful to Professor I. Alexeff for a very thorough and critical discussion of this paper. We would 
like to acknowledge helpful discussions with Professor G. Joyce. In addition we thank Alfred Scheller 
for construction of much of the apparatus. 
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